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ABSTRACT 

We positionally match a sample of infrared-selected young stellar objects (YSOs), iden¬ 
tified by combining the Spitzer GLIMPSE, WISE and Herschel Space Observatory Hi-GAL 
surveys, to the dense clumps identified in the millimetre continuum by the Bolocam Galactic 
Plane Survey in two Galactic lines of sight centred towards I = 30° and / = 40°. We calcu¬ 
late the ratio of infrared luminosity, Air, to the mass of the clump, Mc\ump, in a variety of 
Galactic environments and find it to be somewhat enhanced in spiral arms compared to the 
interarm regions when averaged over kiloparsec scales. We find no compelling evidence that 
these changes are due to the mechanical influence of the spiral arm on the star-formation 
efficiency rather than, e.g., different gradients in the star-formation rate due to patchy or in¬ 
termittent star formation, or local variations that are not averaged out due to small source 
samples. The largest variation in LiR/Mdump is found in individual clump values, which fol¬ 
low a log-normal distribution and have a range of over three orders of magnitude. This spread 
is intrinsic as no dependence of LiR/Mdump with Mdump was found. No difference was 
found in the luminosity distribution of sources in the arm and interarm samples and a strong 
linear correlation was found between Air and Mdump- 
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1 INTRODUCTION 

It is clear from many images of “grand design” spiral galax¬ 
ies that star formation and spiral arm topology are closely inter¬ 
linked. However, the precise relationship between the physics of 
the spiral arms and the nature of their ensuing star formation is 
not yet clear. The influence on the spiral arms could either be 
that they are triggering o ne, or more, o f the main evolutionary 
stages of star formation (iRobertsI 1969|), or are organising fea¬ 
tures for the star-forming material dElmegreen & Elmegreen|[l986l : 
IVogel. Kulkarni. & Scovillelll988h . The main stages of star forma¬ 
tion are represented by: the conversion of atomic gas to molecu¬ 
lar clouds, molecular clouds with internal dense clumps and the 
clumps forming cores from which the single stars form. 

The evidence of the first of these potential influences 
is ambiguous. Increases in the efficiency of cloud formation 
from the neutral interstellar medium have been both observed 


I Hever & Terebevll998h and predicte d ( Dobbs. Bonnell. & Pringld 
I2OO6I) but the results of iFwle et al.l (l201(lh contradict this. They 
find no increase in that efficiency, or the star-formation efficiency 
(SEE) in or out of the spiral arms of two spiral galaxies, NGC 
628 and NGC 5194. There is also no observed difference in the 
fraction of the total molecular gas mass in the form of dense, star¬ 
forming clumps within clouds located in and out of spi ral arms in 
the Milky W ay, when averaged over kiloparsec scales dEden et alJ 
I2OI2L l2013h . There are, however, large cloud-to-cloud variations 
in the dense-clump fraction. There are large increases in the star- 
formation-rate density (Esfr in units of Mq yr“^ kpc“^) by fac¬ 
tors of between 3 and ~ 30 associated with spiral arms in the 
Galaxy. However, no more than 30 per cent of this results from in¬ 
creases in the integrated luminosity of the YSOs per unit molecular 
gas mass, while the rest is due to source crowding. Eurther, much 
of the residual increase in L/M can be ascribed to the inclusion 
of extre me star-forming co mplexes such as W49 and W51 in the 
samples dMoore et alj|2012h . 
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This source crowding suggests that the dominant effect 
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is due to the organising nature of the spiral arms. However, 
IPobbs. Burkert. & Pringle! ( 1201 ih found that the spiral arms al¬ 
low larger molecular clouds to form, by delaying clouds in their 
orbits, which may affect the mass function of cluster s that form 
within them (e.g. IWeidner. Kroupa, & Bonn^ l201(t) . The star- 
formation rate or efficienc y in the higher-colum n density clouds 
are not increased, however (IKrumholz et al ■l201(]h . The spiral arms 
may also be prolonging the lifetime of the clouds, lengthening the 
timescale for star formation to occur, increasing the star-formation 
efficiency (SFE). A study using the Galactic Ring Suryey (GRS; 
Ijackson et alj|2006h found that the clouds i n the interarm regions 
dissipate qu ickly l Roman-Duytdet ai]|2009h with models also pre¬ 
dicting this dPobbs & Pringlell2013 1. 


The aim of this study is to investigate the efficiency with which 
the dense clumps form stars and measure the variation of that effi¬ 
ciency from clump to clump and on kiloparsec scales, looking for 
changes associated with environment. To do this, it is necessary to 
use infrared wavelengths that trace the YSOs embedded within the 
clump. We therefore aim to combine infrared surveys that trace the 
embedded stellar pop ulation of YSOs an d the Bolocam Galactic 
Plane Survey (BGPS: lAguirre et ^1201 ih to calculate the ratio of 
the infrared luminosity to the clump mass, an analogue of the SFE, 
as a function of Galactocentric radius and Galactic environment. 
This stage of the hierarchical, evolutionary star-formation process 
(which converts atomic gas to molecular clouds to dense clumps 
to YSOs) is the last to be investigated. Our previous work has 
found factors of ~ 2-3 increases in the ratio of YSO luminosity to 
molecular-cloud mass, in clouds associated with spiral arms, but no 
enhancement in the average d ense gas mass fraction dMoore et alj 
I 2 OI 2 I : lEden et aDl2012L l2013l) . We would therefore expect to find 
similar factors of ^ 2-3 increases in the average clump to YSO 
conversion efficiency, measured by Lm/Mciump, on similar scales, 
although the results may be sample-dependent in detail. 


The t wo slices of the p l ane s tudied are the same as those 
outlined in lEden et ^ ( 1201 21 l2013h . the I = 30° region spanning 
/ = 28?50-31?50 and the / = 40° region covering Galactic longi¬ 
tudes of / = 37?83-42?50. The 1 = 30° region extends to Galactic 
latitudes of \b\ ^ 1°, whereas the Z = 40° region only covers lati¬ 
tudes of \b\ ^ 0?5. The / = 30° region has three populations within 
it: foreground material which is possibly associated with the Sagit¬ 
tarius arm, the tangent region which is the tangent of the Scutum- 
Centaurus arm and the end of th e Long Bar, which hosts th e ex¬ 
treme star-forming region W43 dNguven Luong et aPbOllf) . and 
the background spiral arms, which are the Sagittarius and Perseus 
arms. The / = 40° region contains two passes of the Sagittarius arm, 
the outer edge of the Scutum-Centaurus tangent and the Perseus 
arm, as well as the interarm regions between the spiral arms. 


A four-armed model of the spiral ar m is assumed, as refined 
by m u ltiple tracers, from H II regions (e.g . iGeorgelin & GeorgelinI 
1 19761: IPaladini. Davies. & De Zottil |2004) to the Hi distribution 
( Levinej_BlitZj&Heif^ 200^ to combinatio ns of tracers (e.g. 
IHou. Han. & Shil l2009l . lUrauhart et alj|2013al l9l. The four spiral 
arms are the Norma, Scutum-Centaurus, Sagittarius and Perseus 


In Section 2, we outline the data sets which are used for this 
project, with Section 3 detailing the reasoning behind the wave¬ 
length choices and the calculation of YSO luminosities, whilst in 
Section 4 we present the results and analysis and Section 5 is a 
discussion of the results. A summary and the concluding remarks 
follow in Section 6. 


2 DATA SETS 
2.1 Infrared data sets 

The Galactic Legacy Infrared Mid-Plane Survey Extraordinaire 
(GLIMPSE) was a Spitzer Space Telescope Legacy Program, 
which (and with subsequent follow-ups) mapped t wo-thirds of the 
Galactic plane as part of thr ee separate surveys dBeniamin et alj 
l2003l:IChurchwell et al.ll200^ . The / = 30° and Z = 40° regions fall 
in the section observed during GLIMPSE I, which covered the lon¬ 
gitudes of |Z| = 10° -65°. The survey mapped this region at wave¬ 
lengths of 3 .6, 4.5, 5.8 and 8 .0 pm using the Infrared Array Cam¬ 
era (IRAC: Ipazio et al.|l2004l) at spatial resolutions of <2". The 
sensitivity limits at each wavelength were 15.5, 15.0, 13.0 and 
13.0 magnitudes, respectively with the high reliability (^ 99.5 per 
cent complet eness) Point Source Ca talogue containing ~ 30 mil¬ 
lion sources dChurchwell et al.ll2009h . _ 

The Wide-field Infrared Survey Explorer (WISE: IWright et aP 
I 2 OIOI 1 mapped the whole sky at 3.4, 4.6, 12 and 22 pm with an¬ 
gular resolutions of 6.1", 6.4", 6.5" and 12", respectively. 5 -ct 
point-source sensitivities of 0.08, 0.11, 1 and 6 mjy were achieved. 
A saturation limit of 330 Jy is found for the brightest sources at 
22 pm, but this is considerably deeper th an MSX at 21 pm, ~ 2.7 Jy 
dPrice et al.ll200ll : lLumsden et alj2013h . 

Hi-GAL (Hers chel infrared Galactic Plane Survey; 
iMolinari et al.ll2010^ f9l. a Herschel Space Obsen’atory Open-time 
Key Project, has mapped the entire Galactic plane at 70, 160, 
250, 350 and 500 pm at spatia l resolutions of 4"-40" using 
the PACS dPoglitsch et alj|2008h and SPIRE dCriffin et alj|2008h 
cameras. The 90 per cent completeness limits at each wavelength 
are 0.5, 4.1, 4.1, 3.2 an d 2.5 Jy at the five wavelengths, respectively 
dMolinari et al.ll20103l . for the Z = 30° region. 

A custom source extraction algorithm , CuTEx (Curvature 
Thresholding Extractor: IMolinari et al.ll^llh . was designed to ob¬ 
tain source catalogues in each of the five Hi-GAL wavebands. 
CuTEx differentiates the image along 4 directions, X, Y and the 
two diagonals twice. Negatives found in the second differential im¬ 
age are taken to be sources as these are the peaks against the back¬ 
ground as long as they are coincident with zeroes in the first differ¬ 
ential image. A source is identified when it is above the threshold 
in all 4 differential directions. The identified sources are then fitted 
with 2D elliptical Gaussian profiles, centred at the positions above 
to extract the integrated fluxes. 

These infrared surveys and their extracted catalogues provide 
almost complete coverage of the mid-infrared to the sub-millimetre 
wavelength range, covering the whole SED of embedded YSOs. 


2.2 Bolocam Galactic Plane Survey 

The BGPS mapped the two regions, Z = 30° and Z = 40°, as part 
of the 133 deg^ northern Galactic Plane section of the survey 
(Z= -10?5to90?5). The BGPS surveyed the continuum at 1.1 mm 
{o = 21lA GHz), with an rms noise of 11-53 mJy beam~^ and an 
effe ctive angula r resolution of 33". 

lEden et ^ ( 1201 21 1201 3l) derived distances, and as a result, 
masses and positions with respect to spi ral arms or with Galacto - 
centric radius to 989 catalogued sources dRosolowskv et alj|201(lll . 
793 and 196 for the Z = 30° region and Z = 40° region respectively. 
The masses found for the BGPS sources were predominantly in 
the rang e of masses associa ted with molecular clumps (10^ ®- 
10® M 0 ; IPunham et al.ll20llh . large, dense, bound regions within 
which stellar clusters and large systems form. 
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The 99 per cent completeness limits in the two fields are 
450 Mq and 550 Mp at 14kpc in the 1 = 30° and / = 40° regions, 
respectively dRosolowskv et alj|201(lh . 


3 SOURCE LUMINOSITY CALCULATIONS 
3.1 Infrared wavelength choice 

A YSO spectral energy distribution (SEP) has two components 
dMolinari et alJl200^ : ICompiegne et al ll2010h . The infrared part of 
the SED comes from the embedded YSO (if already nuclear burn¬ 
ing) and/or the accretion luminosity, while the larger cold envelope 
within which the YSO/HlI region is embedded emits strongly at 
the sub-millimetre and millimetre wavelengths. The latter compo¬ 
nent may be physically independent of the star-formation rate or 
star-forming content of a given clump until the clu mp material be¬ 
gins t o be cleared in the lat er evolutionary stages dMolinari et alj 
l2008l:IUrQuhart et alj|2013d) . 

By fitting an SED to the Hi-GAL data, the suitability of in¬ 
cluding the longer infrared and sub-millimetre wavelengths in a 
measurement of the YSO luminosity can be tested. These SEDs 
were fitted to the band-merged Hi-GAL sources that were position- 
ally matched with the I = 40° BGPS sources. The I = 40° region was 
used as it has the lowest confusion levels of the two fields. 

Before fitting a modified blackbody to the Hi-GAL SEDs, 
we scaled the SPIRE fluxes (Js 250 pm) by the ratio of the beam- 
deconvolved so urce size and beam size at the considered SPIRE 
wavelength (cf. iNguven Luong et^l201 ih . This corrects for the 
increase in the observed source size, an d consequently the f lux, 
at the longest Herschel wavelengths (e.g. iGiannini et alj|201^ . as 
well as for the larger beam size. The sub-millimetre fluxes are likely 
to be dominated by a cold envelope that is physically much larger 
than the embedded YSO. The longer-wavelength sources also tend 
to contain multiple shorter-wavelength detections within the larger 

beam. _ 

_ The fitting strategy is that as described in IGiannini et alj 

( I 2 OI 2 I) . with the modified blackbod y expressions and adapted con¬ 
stants the same as lElia et"^ (l2013h . The free parameters are inte¬ 
grated flux, which is then converted to mass, M, and temperature, 
T, while the dust opacity exponent, p, is kept constant at 2. The 
fits do not include the 70 pm point as the greybody envelope may 
be contaminated by the emission of its protostellar content and this 
point generally lies above the fit. 

The masses derived from the SED fits were compared to those 
calculated from BGPS data, both using a fixed temperature of 14 K 
(see figure 3 of lEden et al.ll201^ and with temperatures obtained 
from the SED fits. The result is shown in Fig.[T] The 1:1 ratio line is 
plotted over the distribution and it can be seen that the mass derived 
from the SED fitting is tracing the same mass as the l.l-mm BGPS 
flux. This distribution also supports the use of a single temperature 
in the previous studies. 

Since the two methods trace the same mass, it is safe to as¬ 
sume that the luminosity traced by the Hi-GAL data at wavelengths 
covered by the SED fit (160-500 pm) is dominated by the clump 
and the cold gas/envelope. This is not useful if trying to identify 
the clump-to-stars stage of the evolutionary sequence and measure 
its efficiency via L/M as M and L are not independent, and if 
they are calculated from the same data, M is essentially L/T. As 
a result, it is more pertinent to use shorter wavelength data, namely 
the 8 pm GLIMPSE data, 12 and 22 pm WISE data and the 70 pm 
Hi-GAL data. These will trace the luminosity produced by the ac¬ 
cretion disk and the embedded YSO, i.e. the stellar component, and 



Figure 1. Comparison of SED derived clump masses and masse s calculated 
from the BGPS 1.1 -mm flux density alone bv lEden ^ ^ j2013h . The black 
crosses represent the BGPS masses calculated using the SED-derived tem¬ 
perature, whilst the red circles use a flat temperature of 14 K. The 1:1 line 
is indicated. 


the calculated luminosity and mass will be relatively independent. 
The wavelengths shorter than 8 pm are not used as they contribute 
very little to the luminosity. 


3.2 Source associations 

For a source to be considered a candidate YSO, it had to be cata¬ 
logued in all 4 wavebands at 70 pm or shorter (i.e. 8, 12, 22 and 
70 pm). The wavebands were combined by positional association. 
Starting at the largest beam size, an association in the waveband 
with the next largest was established if a source was found in that 
band within the beam radius. The order of the association was as 
follows: the WISE data at 6.5” and 12” for the 12 and 22 pm data 
respectively (the WISE all-sky survey contains an already band- 
merged catalogue, so only sources with a flux, and an error on 
the flux, were used), the 70 pm Hi-GAL data with an aperture of 
4” then the Spitzer GLIMPSE 8 pm with the smallest aperture of 
1 . 9 ”. Table Ul gives the numbers of sources at each wavelength. 
The matching process resulted in multiple associations as data with 
smaller beam sizes were introduced. It is usually assumed that only 
one of the smaller-beam-size sources is the “true” counterpart to 
the larger-beam-size components and the determination of the ef¬ 
fect of an incorrect association where there is a choice is discussed 
in Section 3.4. A colour cut was applied that corr esponded to the 
colou r criteria used by the Red MSX Survey (RMS: lLumsden et alj 
l2013h . which was T 22 ^.m > Ts/im- This limits the sample to the 
embedded YSO stage similar to RMS, resulting in an infrared se¬ 
lected sample of 1,641 sources, 1,050 and 591 in the 1 = 30° and 
/ = 40° regions respectively. 

These associations were then positionally matched to the l.l- 
mm BGPS sources, with only those sources that fell within the ex¬ 
tent of a BGPS source counted, as catalogued bv iRosolowskv et ^ 
( l2010h . This allows the distance of the BGPS source to be assigned 
to the infrared (IR) source and only allow for IR sources to be 
counted which are part of a clump. These IR-BGPS associations 
resulted in a total of 241 infrared sources in the I = 30° region and 
107 in the / = 40° region with assigned distances. 

In total, 298 BGPS clumps are associated with YSOs, 29 per 
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Table 1. Source counts at each wavelength in the / = 30° and / = 40° regions. 


Survey 

Wavelength 

/ = 30° 

Sources 

/ = 40° 

Sources 

GLIMPSE 

8 |4.m 

204,367 

163,848 

WISE 

12 & 22 (4.m 

57,506 

47,179 

Hi-GAL 

70 pm 

4,578 

1,469 

BGPS 

1.1 mm 

806 

229 


cent of the total 1,035 BGPS sources. Icsengeri et al.l i2014h found 
a 30 per cent association when matching ATLASGAL and WISE 
sources. Of these 298 clumps, the vast majority are associated with 
a single infrared source (250 or 84 per cent), with only two clumps 
found to be associated with 3 infrared sources, one in each of the 
two regions. 

To estimate the number of chance alignments, a Monte Carlo 
simulation was conducted, in which 1,641 sources were distributed 
over the / = 30° and / = 40° regions. The random distribution of sim¬ 
ulated sources accounted for the true clustering scales of the IR 
sources by maintaining the relative position of sources found within 
a BGPS beam of each other. These sources were then matched with 
the BGPS source catalogue. 100 repeats found that the percent¬ 
age of sources aligned with BGPS detections had a mean value of 
0.8 ±0.3 per cent, compared to the actual association alignments 
of 14.7 per cent. 

3.3 Galactocentric radius distributions 

Using the assigned velocities and distances from lEden et alJ ( l2012L 
l2013f l for the BGPS clumps, and hence the YSOs, the distribution 
of the clumps and YSOs as a function of Galactocentric radius can 
be produced. These are displayed in the upper panel of Pig.|2] The 
ratio of these distributions is presented in the lower panel. This 
quantity gives the number of YSOs per dense clump and can be 
used as an analogue of the star-formation efficiency. The absolute 
luminosity-to-mass ratio as a function of Galactocentric radius is 
presented in Pig.|7]and discussed in Section 4.3. The peaks in both 
distributions occur at the positions of the Scutum-Centaurus tan¬ 
gent and W43. This is not surprising considering that the line of 
sight is down the spiral arm and W43 is the most prominent star¬ 
forming region in either of the lines of sight. 

The number count ratio is relatively flat across most Galacto¬ 
centric radii, with almost all the values for each 0.5 kpc bin falling 
within 1- or 2-a of the mean. The mean value is found to be 
0.46 with a standard deviation of 0.20 and a standard error on the 
mean of 0.07. The exceptions, however, are the bins at Galacto¬ 
centric radii of 4-5 kpc, with values ~ 5-a below the mean. These 
radii are associated with the S cutum-Centaurus tangent as well 
as the W43 star-forming re gion. iMotte, Schilke, & Lda l l2003h and 
iNguven Luong et^ ( 1201 ih interpret this as evidence that the star- 
formation rate (SFR) will be higher in the future in W43 due to an 
excess of starless clumps. 

3.4 Luminosity calculations and the effect of incorrect 
association 

Two methods were used to calculate the YSO luminosities: an inte¬ 
gration usingjhebapeziummle method and the online SED fitting 
tool of iRobitaille et alj l l2007tl . The latter makes use of a library of 
200,000 models of YSOs. A detailed description of the models can 
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Figure 2. Upper panel: Galactocentric radius distribution of the BGPS 
clumps (grey filled histogram) overlaid with the YSOs (blue dotted his¬ 
togram). Lower panel: Ratio of the two distributions, with the mean value 
of 0.46 indicated by the dashed line. The errors are derived using Poissonian 
statistics. 

be found in iRobitaille et al.l ( l2006h and iRobitaille et al] ( l2007h . The 
trapezium integration assumes no physical model but may underes¬ 
timate the total luminosity as there is no correction for the missing 
wavebands or contributions to the IR flux from the cold envelope 
component. The SED fits are based on models that may not be ap¬ 
propriate but provide the above extrapolation based on a spectral 
shape that is generally plausible. 

The distribution of the Robitaille SED calculated luminosities 
compared to the trapezium rule luminosities is shown in Eig. 
The comparison shows the 50 per cent limits of the distribution of 
L value ratios and 63 per cent of sources are found within these 
limits. The distribution is skewed to excess luminosity in the model 
SEDs. As can be seen from the right-hand panel of Eig. [4] SEDs 
fitted to these sources tend to have a significant component of emis¬ 
sion at near-IR, or sub-millimetre, wavelengths, where we have no 
data constraints. This is symptomatic of the sources with a SED 
fit which overestimates the IR luminosity, a large component from 
either the sub-millimetre or near-IR. The sources with a greater 
trapezium integration are over-luminous at 70 pm. This makes the 
Robitaille fits unsuitable for those purposes as we are attempting 
to determine the IR luminosity independent of the sub-millimetre 
clump component. Fig. |4] shows an example SED fit with a consis- 


© 2015 RAS, MNRAS POQ.fTHTTI 



















































SFE in spiral arms 5 



log,oBolometric Flu* from Robitoille Fitter/L© 


Figure 3. The comparison of IR-BGPS sources showing the correlation be¬ 
tween the trapezium rule integration and the Robitaille model luminosity. 
The 1:1 line is shown as the dashed line whilst the 50 per cent limits are the 
blue dot-dash lines. 


tent trapezium-rule luminosity and one with a significantly overes¬ 
timated Robitaille luminosity. 

As a result, the luminosities of the YSOs used for this study 
were estimated via the trapezium rule. This however introduced a 
potential systematic error to the calculations since any emission at 
wavelengths < 8 p.m or > 70 [xm is not accounted for. However, the 
absolute values of the luminosity are less important than the com¬ 
parative values_and_si^ biases are assumed to be similar across the 
sample. FVeneziani et al.l l l2013ll estimated bolometric luminosities, 
rather than infrared luminosity, via a similar method to the trapez¬ 
ium rule, using those luminosities to produce star-formation rates 
and L — Mdump relations for the two Hi-GAL science demonstra¬ 
tion fields. 

As mentioned above, multiple GLIMPSE and Hi-GAL 
sources were often found to be associated with a single WISE 
source, or vice versa, due to the different spatial resolutions. We 
therefore investigated the effect of making an incorrect choice of 
multiple components on the final calculated luminosity. 

By fitting SEDs to all the potential infrared source matches 
for a given BGPS source in a subset of the sample, we investigated 
the uncertainty in luminosity values from the choice of source. Al¬ 
though some match combinations produced drastically better fits 
than others, the resulting luminosities were found to vary by less 
than 10 per cent. It is therefore clear that, for these purposes, there 
is no need to be overly concerned with identifying the correct coun¬ 
terpart and that the luminosity is dominated by wavebands (22 pm 
and 70 pm) in which there are few or no choices. Use of the trapez¬ 
ium rule produced similar results, with luminosities from each po¬ 
tential set of matches that were consistent within the measurement 
uncertainties. However, the absolute luminosities produced by the 
two methods were different. 


4 RESULTS AND ANALYSIS 
4.1 Star-formation efficiencies 

The star-formation efficiency, or SFE, considered here, one of sev¬ 
eral analogous variants, is the fraction of dense, dust-continuum 
traced clump material that has been converted into stars over some 


timescale. As with the clump-formation efficiency (CEE; the rati o 
of clump mass to cloud mass) analysis in lEden et alj l l2012Ll2013h . 
the SFE within the separate populations of clumps along the two 
lines of sight can be estimated. The SFE is the ratio of the luminos¬ 
ity produced by embedded young stellar objects to the mass of the 
clumps. 

The ratio of infrared luminosity of YSOs to clump mass, is a 
measure of: 


L.tar _ 1 p ^ 

-^^clump -^^clump J Q dt 


where dL/dt is the instantaneous star-formation rate in terms of L. 
A high value for L/M can indicate either a high star-formation 
rate or a longer timescale. Equating L/M with the SFE there¬ 
fore requires a few assumptions. For precise linear dependence, 
dL/dt must be proportional to dM/dt which requires that the 
stellar IMF is fully sampled in all star-forming regions, up to the 
maximum stellar masses, i.e. “sorted sampling” filling the IMF 
jWeidner & Krout^l2006h . If, as se ems likely, this is not the case 
and the IMF is filled stochastically ( lElmegreenll2006l) . then L/M 
still has a direct dependence on the SFE but it may be non-linear. 
Then, an increase in SFE must produce an increase in L/M but 
an observed increase in L/M can also be the result of the for¬ 
mation of a larger cluster at the same efficiency (L oc for 
clusters). Such an increase in L/M cannot be accounted for ob- 
servationally until it is possible to count all the individual stars 
within a cluster. Another source of non-linearity may be intro¬ 
duced if the SFE is high, in which case, Mdump is no longer con¬ 
stant, bu^_SFE_va^es_appea£_to_be_generalI^ess than 30 per cent 
(e.g. iLada & Ladall2003l : [Johnston et alj|200^ and so Mdump is 
taken to be constant over the lifetime of a star-forming event as 
traced by mid- and far-infrared detections. In principle, L/M also 
increases with time and the SFE is usually defined with respect 
to some characteristic timescale; however, the stage of massive 
star formation trace d by mid-infrared emission is short (less than 
a few X 10® years: lOavies et al.ll201 ll ; iMottram et alj|201 lah and 
our L/M determination can be taken as a snapshot of the current 
SFE. Vari ations in L/M can also be used as potenti al evolution 
indicators dMolinari et alJl200^ : IUrauhart et alj|2014bl) although L 
and Mdump tend to show a linear relationship with little evidence 
of evo lution in L/M during the IR-bright phase dUrquhart et alj 
l2014hh . 


The clump populations can be divided into those belonging 
to the foreground spiral arms, the tangent region consisting of the 
Scutum-Centaurus tangent and the end of the long bar, and the 
background Sagittarius and Perseus arms in the / = 30° region and 
the spiral arm and interarm regions in / = 40°. The 1 = 30° region 
populations were separated using the distribution of the distances to 
th e molecular cloud s fromiRom an-Duv al et alj d2009l) . as outlined 
in lEden et al.l d2012h . lEden et al. (2013) details how the molecular 
clouds, and hence, the BGPS clumps were assigned to spiral arm 
and interarm pop ulations in th e 1 = 40° region using the Galactic 
rotation model of IValleel dl995l) . 

Using t he masses derived for the BGPS sources bv lEden et alj 
d2012Ll2013l) and infrared luminosities for the associated YSOs ob¬ 
tained using the trapezium-rule integrations, we calculate the inte¬ 
grated luminosity and clump mass in each of the populations in the 
1 = 30° and / = 40° fields and use the ratio, L/M as an estimate of 
the SFE averaged over the scale of the field. The values for each 
region are displayed in Table[2 BGPS-traced clumps are included 
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Figure 4. Example Robitaille SEDs with the trapezium-rule integration limits mai'ked by the two blue lines. The SED-fit luminosity is the total area under 
the curve. The black points ai'e the data points to calculate the best fit, the black line represents the best fit with the next best 9 fits represented in grey. The 
dashed line shows the best model if a stellar object is assumed with no circumstellar extinction. Left panel: A source which exhibits similar luminosities for 
the Robitaille and trapezium-rule fits is shown with a trapezium-rule luminosity of 1.4 x 10"^ Lq and Robitaille-fit luminosity of 2.4 x 10^ Lq. Right panel: 
A source which exhibits an extreme difference in luminosity between the Robitaille and trapezium-rule fits is shown with a trapezium-rule luminosity of 
S.lxlO^Lo and Robitaille-fit luminosity of 4.0 X 10® Lq. 


whether or not they contain IR sources since they contribute to the 
total L/M. 

Th e results for t he / = 30° region differ from the results that are 
found in lEden et alj ( 12012h who found consistent clump formation 
efficiencies in the tangent and background region and an increase in 
the foreground. The SFEs for each region depart from each other by 
factors of ~ 2, that are significant within errors. The / = 40° region 
shows that arm population has a value twice that of the interarm. 
The implications of this are discussed in Section 5.1. 

The uncertainties on the S FEs are a combination of the uncer¬ 
tainties in BGPS flux densities dRosolowskv et al]|2010h . the vari¬ 
ance of the temperature distribution of the source populations, the 
catalogued flux uncertainties on the GLIM PSE and WISE sources 
dChurchwell et alj2009l:rWright et alj20T^ . as well as the recently 
released Hi-GAL catalogues (Molinari et al., in prep) containing 
the source flux uncertainties. 

In summary, the average SFEs for the spiral arms are enhanced 
compared to the interarm regions by a factor of ~ 2. However, there 
is a suppression towards the W43 star-forming region back to the 
interarm level. The other spiral arms, especially the Perseus arm, 
show an elevated SEE. 


4.2 Luminosity functions 

The luminosities of the YSOs were used to plot the luminosity 
functions (LFs) of the Scutum-Centaurus tangent region and of 
the other spiral arm regions. The spiral-arm subsample includes all 
YSOs found in the background region of 1 = 30° and the Perseus 
and Sagittarius YSOs from / = 40° as the background region con¬ 
sists of Sagittarius and Perseus YSOs. All sources with a helio¬ 
centric distance of less than 2 kpc are omitted, in order to remove 
local sources that may contaminate the results and avoid including 
a large number of low-luminosity YSOs below the completeness 
limit at Galactocentric radii of 8 kpc. No interarm LF is made due 
to the lack of sources. 

The LFs of these two populations are presented in Fig. [5] 
The plotted quantity, AN/AL, is the number of sources per unit 


luminosity interval in each luminosity bin, with the luminosity 
point representing the median value of the bin. A fixed number of 
sources per bin was used, as opposed to fixed bin widths, in order 
to equalise weights determined from Poisson errors. 

If it is assumed that the data in both populations are distributed 
as a power law of the form AN/AL oc L“, a least-squares fit to 
each LF gives indices of a = — 1.55 ±0.05 and a = — 1.57 ±0.07 
for the tangent and arm regions respectively. Both LFs begin to 
flatten at luminosities ~ log{L /'Lq) = 2.5 which we assume to 
indicate the sample completeness limit. The index values for these 
populations are consistent with each other, which shows that, in 
terms of the resulting distribution of YSO luminosities, the star- 
formation process at the end of the bar is no different to that in 
other spiral arms, despite the inclusion of the W43 region. 

These slopes are consistent with that of the RMS survey 
dLumsden et al.l l2013h. Plotting the lum i nosities of both YSOs 
and H II regions iMottram et al.l 1201 Ibl : lUrouhart et al.l l2014al) 
with equal population bins give indices of a = — 1.48 ±0.04 and 
Q = — 1.60 ± 0.04 for the YSO and H II populations respectively. A 
combined sample had an index of a = — 1.50 ± 0.02. 

Further evidence for the similarity of the star formation in the 
two populations is that a Kolomogorov-Smirnov (K-S) test shows 
a 27 per cent probability that they are drawn from the same popu¬ 
lation, which is only a 1-cr departure. Thus we are unable to reject 
the null hypothesis that they are drawn from the same population. 
Even though no LF was produced for the interarm regions, by per¬ 
forming a K-S test with the interarm sources, it showed a 23 per 
cent probability that they are drawn from the same population so 
there is no evidence for a variation in the current data. The cumu¬ 
lative distributions of the three samples are displayed in Fig.|^ and 
the similarity of the three distributions reflects the K-S test results. 
It is worth noting that the interarm sample is very small and not a 
strong constraint on the nature of interarm star formation. 
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Table 2. SFEs calculated from the YSO luminosity to clump mass ratio for the populations within the / = 30° and I = 40° regions. 


Region 

Population 

Tir 

(Le) 

No. of 

YSOs 

-^clump 

(M©) 

No. of 
Clumps 

SFE 

(Lq/M0) 


Foreground 

6750 ± 193 

28 

10300 ± 648 

72 

0.65 ± 0.05 

II 

o 

o 

Tangent 

523000 ± 5090 

160 

544000 ± 9990 

615 

0.96 ± 0.02 


Background 

492000 ± 6270 

53 

144000 ± 6880 

106 

2.47 ± 0.09 

II 

o 

o 

Total Arm 

191000 ±2570 

77 

137000 ±9180 

135 

1.40 ±0.10 


Total Interarm 

59700 ± 2880 

30 

83300 ± 7070 

61 

0.72 ± 0.07 
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Figure 5. LFs for the Scutum-Centaurus tangent region and all other spi¬ 
ral arm regions. The dashed lines are fits to the slope of the LF. The eiTor 
bars correspond to \/N counting statistics. Top panel: Scutum-Centaurus 
region. Bottom panel: Spiral arm regions with sources of heliocentric dis¬ 
tances less than 2 kpc omitted from both the I = 30° and I = 40° regions. 

4.3 LjM distribution with Galactocentric radius 

The distribution of L/M as a function of Galactocentric radius is 
displayed in Fig. |7] All sources with heliocentric distances under 
2 kpc are again omitted. A significant peak is found at 6 - 6.5 kpc, 
with an increase by factors of 2.2 and 3.2 over the adjacent bins, 
which correspond to 7.8 -(t and 9.3 -ct peaks, respectively. A less 
pronounced peak is found at 7 -7.5 kpc, with 4.1-a and 3.1-a de¬ 
partures found compared to adjacent bins, which are increases in 
factors of 2.5 and 1.9, respectively. The first peak corresponds to 


Figure 6. The cumulative distribution function of the spiral ai'm, interai'm 
and tangent region samples. The spiral arm sample is represented by the 
black dotted line, the interai'm sample by the blue dashed line and the tan¬ 
gent by the red dot-dash line. The spiral arm sample again has sources of 
heliocentric distances less than 2 kpc omitted. 


the Sagittarius spiral arm, whilst the second is the Perseus spiral 
arm. _ 

iMoore et alj ( 12012h found peaks associated with the Sagittar¬ 
ius and Perseus arms in the L/Mdoud distribution, using the CO- 
traced molecular cloud mass. These peaks were found to be ex¬ 
treme in their study due to the presence of the star-forming regions 
W51 and W49A, respectively. Flowever, when they removed these 
star-forming regions, the Perseus arm had a SFE similar to the lev¬ 
els in the adjacent bins whereas the removal of W51 from the Sagit¬ 
tarius arm reduced the peak relative to the neighbouring bins but the 
peak remained. 

This study has neither the W49 or W51 star-forming regions 
within the two lines of sight investigated . As a result, the di stribu- 
tion as seen in Fig.|7]should reflect that of iMoore et alj ( l2012h . once 
those regions are omitted as the Md umu/Mdoud ratio is co nstant 
(~ 5 - 10 per cent) over large scales dEden et al.ll20lll2013ll . The 
peak associated with the Sagittarius spiral arm is still elevated com¬ 
pared to the adjacent bins, whereas the Perseus arm bin at 7 kpc has 
a marginal peak comp ared to the bins either side, hence reflecting 
the lMoore et al.l ( 12012h distribution. 

No peak is found at 4 kpc, the radius associated with the 
Scutum-Centaurus tangent and the extreme-star-forming region, 
W43, suggesting that the current star formation in W43 is unspec¬ 
tacular in terms of the YSO luminosity being produced per unit gas 
mass. 
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Figure 7. The total L/M as a function of Galactocentric radius within the 
studied regions centred at / = 30° and I = 40°. The Galactocentric radii as¬ 
sociated with spiral arms are labelled. Sources with heliocentric distances 
less than 2 kpc are omitted, to remove local sources at Galactocentric radii 
of ~ 8 kpc. 

4.4 L/M ratios of individual clumps 

The L/M ratio of individual BGPS clumps can be calcu¬ 
lated and the distribution of these values ranges from 0 to 
14 OLq/M 0. The mean value is 5.24 ± 0.70 Lq/Mq with a median 
of 1.72 ±1.14 Lq/Mq. 

The distributions of L/M values in the separate regions listed 
in Table [2 are statistically indistinguishable and can be assumed 
to be a single population. We can therefore combine these samples 
and the distribution of the total L/M ratios for individual clumps is 
shown in Fig. [8] The total distribution, as with the clump-formation 
efficiencies i mEden et al.l(l2012Ll2013l) . can be considered a lognor¬ 
mal distribution as Anderson-Darling and Shapiro-Wilk normality 
tests gave probabilities of 0.434 and 0.563, respectively, that the 
distribution is lognormal. The presence of a log-normal distribution 
indicates that the variation in L/M is likely to be due to the random 
variations in a collection of individual processes that multiply to¬ 
gether to produce the observed values, and this suggests that those 
values in the wings of the distribution, both with low and high L/M 
are not the result of fundamentally different processes. No correla¬ 
tion of L/M with clump mass was found. Therefore the spread of 
values in Fig.[^is intrinsic and the result of natural variations in the 
star-formation process within clumps that are independent of this 
basic property. 

The lognormal distribution can be fitted with a Gaussian that 
has a mean of 0.08 and a standard deviation of 0.57 dex. With these 
parameters, three clumps can be found to have IR luminosity-clump 
mass ratios that are greater than the mean by more than 3 -ct. With 
the probabilities of these events, only one source at these signif¬ 
icances would be expected in a sample of this size. These three 
clumps are all found in the Scutum-Centaurus tangent in the / = 30° 
region. These clumps have BGPS IDs of 4497, 4561 and 4617 with 
masses of 104, 197 and 97 Mq, respectively. These all have one as¬ 
sociated YSO which have luminosities of log{L/L q) = 3.94, 4.15 
and 4.13. 

The absolute SEE values (Mstar/ATdump) can be calcu¬ 
lated for all the clumps which have an associated luminosity of 
loq(L/hpi) > 3.01 by using the stellar luminosities outlined in 
iMottram et alj ( 1201 lah for spectral types B3 and greater, calculated 


50 



- 2-10 1 2 
log,|,(Lummosity/Clump Moss (L0M0"')) 


Figure 8. The ratio of YSO luminosity to BGPS clump mass for all clumps 
with at least one associated YSO in the 1 = 30° and 1 = 40° regions. The bin 
size is 0.2 dex. 

by using a 1 Myr isochrone to local OB stellar popula tions using 
the stellar rotation models of iMevnet & Maed^ J2003h . The range 
of these SFEs is from 0.1 to 14.3 per cent with a mean value of 
1.9 ±0. 3 per cent an d a median of 0.9 per cent, consistent with 
those of ISmithI (l2014ll . As these are most likely to be the sites of 
cluster formation, the SFEs calculated here are lower limits. This is 
only the star-formation efficiency derived from the brightest source 
and likely missing most of the lower-mass YSOs. 

As we have determined the L/M ratios of the individ¬ 
ual clumps that contain a YSO, the range can be compared to 
those found for molecular clouds and whole spiral arms, to de¬ 
termine the scale at which the most drastic variations are occur¬ 
ring. The range of L/Mdump for the individual clumps that have a 
YSO is 0.04±0.01-140± 22LqMq. The range of L/Mdoud 
is O.O2±O.Ol-68±13Lo/M0 and L/M^,m is 1.00±0.02- 
2.9±O.2L0/M0. 


5 DISCUSSION 

5.1 Star formation on kiloparsec scales 

lEden et aP (12012112013) found that the dense-clump mass fraction 
in clouds (Mdump/ATdoud) varies by more than three orders of 
magnitude from cloud to cloud (a = 0.57) on spatial scales of 10- 
100 pc. However, on kpc scales and between arm and interarm com¬ 
ponents, Mdump/ATdoud is consistent with being constant, within 
the current uncertaint ies, and is se e mingl y unaffected by the pres¬ 
ence of a spiral arm. iMoore et al.l ( 1201 2h showed that L/Mdoud, 
averaged over kpc scales, increases a little in some spiral arms, 
relative to the nearby interarm regions, but not in others. In par¬ 
ticular, the Scutum-Centaurus arm tangent shows no increase in 
L/Mdoud, despite being the location of a large enhancement in the 
surface density of the star-formation rate (Esfr). Most of the large 
increases in Esfr associated with spiral arms are therefore due 
to crowding, rather than a change in the star-forming capacity of 
clouds. The residual enhancements in L/Mdoud that were found in 
other arms are largely due to the inclusion of extreme star-forming 
regions such as W49A and W51 in the samples examined. 

From the current work (Fig. |7]l, we see that there are also 
increases in L/Mdump associated at least with the Sagittarius 
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Figure 9. The L/M as a function of Galactocentric radius within the stud¬ 
ied regions centred at / = 30° and / = 40° with the non-star-forming clumps 
removed. Sources with heliocentric distances less than 2 kpc are omitted, to 
remove local sources at Galactocentric radii of ~ 8 kpc. 


arm, and possibly also with the Perseus arm, by factors of about 
2 , compared to the interarm regions, but, again, no enhancement 
in the Scutum-Centaurus tangent. Since we have evidence that 
Afciump/Afcioud is constant on these scales, we expect that any 
variations in L/Mdump in the pr esent study would b e consistent 
with those found in L/Mdoud bv iMoore et alj l l2012h and, allow¬ 
ing for the differences in the samples used in the two studies (the 
current study does not include W49A, for example), this is indeed 
what we find. 

The enhancements by a factor of ~ 2 in the average L/Mdump 
found in some spiral arms (Sagittarius and Perseus) but not in 
others (Scutum-Centaurus) might be explained by different time 
gradients in the star-formation rate in different places. If we re¬ 
move the apparently non-star-forming clumps from the calculation 
of L/Mdump (i-e- those clumps without an IR source), we find 
that a ~ 30 per cent enhancement appears at the radius of the 
Scutum-Centaurus tangent (Fig. with respect to the adjacent 
radii. L/Mdump here is larger by a factor of ~ 3 than the value 
obtained when the starless clumps are included in Mdump- Else¬ 
where within a radius of ~ 6.5 kpc, L/Mdump increases by a factor 
closer to 2 than when starless clumps are not counted. This is con¬ 
sistent with the idea that the Scutum-Centaurus tangent region has 
a larger-than-average reserv oir of dense clumps yet to form stars. 
iNguven Luong et^ ( 1201 ll) predict that the SFR in clouds associ¬ 
ated with the W43 star-forming region, which is within the tangent 
and included in our sample, will be five times greater in the future 
than the current rate. Beyond 6.5 kpc, and including the Perseus 
arm, the value of L/Mdump is less dependent of the inclusion of 
starless clumps, indicating that there are relatively few of these. 

In this scenario, the current gradient of the star-formation rate 
is relatively flat at radii less than 6.5 kpc, including the Sagittarius 
arm. In the Scutum-Centaurus tangent, the star-formation rate is set 
to increase with time and around the Perseus arm, star formation is 
currently relatively advanced and the production rate of clouds and 
of clumps within clouds is not keeping pace with the star forma¬ 
tion. This picture is consistent with the kno wn patchy nature of star 
formation in spiral arms and elsewhere (e.g. lUrouhart et al]l2014bl) 
which naturally suggests that it is also intermittent in time. It may 
also require that the timescale for the infrared-traced star forma¬ 


tion is longer than that for the production of dense clumps, so that 
clump production does not keep up, or that the dense clumps are 
aboriginal cloud features that are not replaced. It is usually assumed 
that since the timescales associated with the star-formation stages 
trace d by infrared and sub-mm emission are a fe w x 10® years or 
less dMottram et ali[20lT^ : lGinsburg et al.Ll2012h . then these trac¬ 
ers give a snapshot of current star formation, with little or no look- 
back time. _ _ _ 

As shown in lEden et aP ( 1201 2l . l2013l) and in Fig. [S] the 
distributions of the clump-formation and star-formation efficien¬ 
cies of different populations of individual clouds all follow sta¬ 
tistically indistinguishable log-normal distributions. Since spiral 
arms generally have more star formation in them, the spiral-arm 
populations are larger and sample further into the wings of the 
L/Mdump probability distribution, and so include more very high 
(and very low) SFE clouds and clumps but, since L{M) for mas- 
sive, star-forming clu mps is linear rather than ZAMS-like (e.g. 
lUrguhart et alj|^14bl) . the sample average of L/Mdump does not 
depend on sample size and should not vary between arm and inter¬ 
arm populations for that reason. This does not rule out the possibil¬ 
ity of some local variance p roduced by the inclu sion of occasional 
extreme objects, as found in lMoore et al.l ( l2012h . 

A physical explanation as to why L/Mdump should be gen¬ 
erally enhanced in spiral-arm clumps is elusive. Star-formation 
timescales may be longer in long er-lasting and generally more mas¬ 
sive clouds dPobbs et al.ll201 ih . but there is no independent evi¬ 
dence to support this. Positive feedback (i.e. triggering) in crowded 
regions might be invol ved, although mo dels suggest that this is 
not a significant effect dPale et al. l2014h and there is no indica¬ 
tion that clumps in spiral-arm clouds are physically different from 
those elsewhere. While a more definitive resolution of this question 
requires studies of larger samples, results so far contain no strong 
evidence of significant variations in SFE on kpc scales that cannot 
be ascribed to the inclusion of extreme objects in the samples and 
hence no indication that large-scale mechanisms in general, and 
spiral arms in particular, are important in regulating Galactic star 
formation. The fact that variations in conditions within individual 
clouds are much more significant, as outlined in Section 4.4, impli¬ 
cates cloud-scale mechanisms such as the initial conditions of cloud 
formation, cloud collisions or local feedback as the key factors. 


5.2 Potential sources of bias 

5.2.1 Distance dependence of L/M and completeness 

One source of potential bias relating to the L/M ratio is the dis¬ 
tance. There are three ways in which distance effects could alter 
the L/M ratio: the nearby clumps have elevated L/M ratios due 
to a greater detection percentage of YSOs, the most distant clouds 
having a deficit due to the YSO completeness dropping off before 
the BGPS completeness or the inverse of this and the lack of BGPS 
completeness pushing the L/M ratio up. 

The detection limits of the YSO-tracing surveys could simul¬ 
taneously increase the L/M ratios nearby and decrease it at further 
distances. A larger percentage of infrared sources would be ob¬ 
served in nearby clumps and clouds as the lower luminosity sources 
would have been detected. This effect would work in reverse at 
larger distances. This effect is not seen in Fig.[T0]which shows a flat 
distribution in the ratio of numbers of YSOs to clumps. The flatness 
of the distribution indicates that the surveys trace similar distances, 
with no enhancement or suppression of the ratio at the furthest dis¬ 
tances. It is worth noting that any YSOs detected in nearby clumps. 
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Figure 10. Ratio of the number of YSOs to BGPS clumps as a function of 
heliocentric distance, with the mean value of 0.49 indicated by the dashed 
line. The en'ors are derived using Poissonian statistics. 


5.2.3 Assumed spiral arm model 

For t his work, and those of the clump-formation efficiency anal¬ 
ysis teden et ^120121 l2013h . a four-armed model of the Galaxy 
is adopted. This model is preferred over contrasting two-arm mod¬ 
els, determined fro m the kinematics of Two Mi cron All-Sky Survey 
(2MASS) sources jFrancis & Anders^l2012h . or t hree-arm mod¬ 
els fr om the distrihutions of giant molecular clouds ISolomon et all 
Il987l) . 

The four-armed model is preferred due to the vast number 
of different tracers whi ch outline four arms, from Hll regions 
(e.g IPaladini et alj[2004l) . to Hi distribution s iLevine et ahlEbodl 
to combinations of tracers (e.g. IHou et a D 1 20091) . However, the 
YSO population of the Galaxy is consistent with a four-arm model 
lUrouhart et alj2014al) . and as we are tracking the YSO population, 
this is the best choice. 

The choice of this model over others would have a bearing on 
the arm/interarm separation of the sources. It, however, would not 
have a bearing on the Galactocentric radius distribution. As such 
it would be clear that on the scale of individual clumps we would 
still see the large variation as seen in Fig. but the variations on 
kiloparsec scales would be averaged out. 



Distonce/kpc 


Figure 11. L/M ratio for individual clumps as a function of distance in the 
1 = 30° and 1 = 40° regions. 


but missed at far distances, will be low luminosity and will not add 
much to T/Mciump- In the reverse case, undetected clumps at fur¬ 
ther distances will not be associated with any bright YSOs. 

The respective completeness limits of the two surveys could 
also alter the resu lt. The BGPS has a co mpleteness limit of ~ 
550 Mq at 14kpc iRosolowskv et al.|[2010h . whilst the YSO lumi¬ 
nosity is limited by the least sensitive survey, Hi -GAL 70 trm, bu t 
has a completeness limit at ~ 675 L© at 14 kpc l lElia et al.ll2(^ . 
These surveys are similarly sensitive out to the furthest distances 
sampled in this study, and the lack of distance bias is shown in 
Fig.E] with a flat distribution found. 


5.2.2 Area of Galaxy studied 

The regions sampled in this study are relatively narrow, in terms 
of the entire Galactic plane. However, it is unlikely that this intro- 
du ces a significant bia s as our results are fully consistent with that 
of iMoore et al.l ( 1201 2h who found similar increases in L/Mdoud 
over a much greater fraction of the Galaxy. 


5.3 Clump mass-luminosity relation 


The relationship between clump mass and YSO luminosity is dis¬ 
played in Fig. There appears to be a correlation, but both quanti¬ 
ties have a strong distance dependence. Therefore, to test the corre¬ 
lation, a partial Spearman correlation test, which takes into account 
the dependence on a common variable (distance in this case), was 
used. This gave a correlation coefficient of 0.56 with a p-value of ~ 
X 10“^ that this correlation occurred by chance for the sample of 
298 BGPS clumps. By selecting a distance-limited sample between 
4 and 6 kpc, a correlation coefficient of 0.52 with a p-value of ~ x 
10 ~® is found, therefore there is no evidence for any significant 
systematic distance effect on the observed relationship. 

Interpretation of the L — Mdump relation is difficult if the L 
and M ciump are derived from the same data (e.g. IVeneziani et all 
l2013h . which is particularly the case if L = Lboi, hence the use 
of Lir instead of Lboi to attempt to remove the dependence of 
both quantities on the sub-mm component of the SED. The linear 
log-log fit of Lir with Lboi in Fig. [3] is the result of YSO SEDs 
being essentially greybody spectra and that the sources are selected 
in the same way i.e. by waveband, with small variations, such as 
the temperature differences reflected in the differing colours of the 
sources and the independence of Lir and Lboi causing the scatter. 

The slope of the linear log-log fit found for thi s distribution is 
0.99 ± 0.04, which is slightly shallow er than that ofiMolinari et all 
( l2008h . but is consistent with that of lUrguhart et al. (l2014bl) who 
found slopes of 1.27 ±0.08 and 1.09 ±0.13, respectively. The 
slope of the linear fit to the distance-selected sample mentioned 
above is 0.89 ± 0.18. 

Eurther to this correlation, the L — Mdump relation can be 
used to derive the ev olutionary state of th e YSO as it approache s 
the main sequence |Molinari et alj l2008l : iGiannetti et all l2013l) . 
iMolinari et alj ( l2008l) highlight two phases in this evolution in 
the L — Mciump plane: the main accretion and envelope clean- 
up phases, using the accelerating accretion rates as postulated by 
dMcKee & Tanll2003h . Once the clump is formed, the gradient of 
the locus in L — Mdump space will begin at Lir = 0, as there 
would be no YSOs present, and so Lir will be zero. The gradient 
will be 1 if L = Lboi and the temperature is constant with L, and 
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Figure 12. Clump mass-luminosity relationship for the matched YSOs and 
BGPS clumps. The I = 30° populations of the foreground, tangent and back¬ 
ground regions are represented by purple plus symbols, black stars and red 
squares, respectively, with the I = 40° arm and interarm sources represented 
by blue triangles and green crosses, respectively. The lower, middle and up¬ 
per black dashed lines represent the = 0.1, 1 and 10 Lq/Mq, 

respectively. The red dot-dash line represents the fit to the data. 

M. As YSOs form, the luminosity will be made up of two compo¬ 
nents, the accretion luminosity and the luminosity of the source. If 
the accretion luminosity, Lacc, dominates the YSO luminosity, then 
this implies a high accretion rate and/or a very early stage prior to 
fusion. The luminosity thus increases at near-constant mass on a 
vertical track in L — Mdump space. Once the accretion finishes, the 
envelope begins to be dispersed. The YSOs find themselves clus¬ 
tering at T/Mciump = 1 Lq/Mq as sources spend a relatively long 
time in the region where accretion is ending and envelope clearance 
has not yet commenced. Both the accretion rate and envelope clean¬ 
up rate increase with time with the accretion rate inc r easing as the 
mass of the protostar increases dMcKee & Tanll2003l : IPavies et alj 
l201lh . The clean-up phase starts slowly as it takes some time for the 
embedded star to disrupt the clump from which it formed. The final 
luminosity is dependent on the initial clump mass, and as a result 
a linear log-log distribution is found, especia lly as the accretion lu - 
minosity is proportional to the stellar mass dMcKee & Tanll2003h . 
This dependence, LocM* ®®, implies the accretion luminosity is 
proportional to the initial clump m ass due to the relatio nship be¬ 
tween stellar mass and clump mass dMcKee & TanllrOO^ . 

Once the clearing is finished, the YSOs move to the left in 
L —Mciump space, which is a very quick phase, hence why no very- 
high-luminosity, low-mass sources are found. By selecting sources 
in the mid-IR, we may be filtering out the envelope-clearing and 
pre-stellar sources. This begins the evolution from the slope of 
unity towards the mass-luminosity relation of either main-sequence 
stars (L oc M®'®) or clusters (L oc M^). 

The distribution also implies that there is a constant 
L/Mciump ratio acting along the spectrum of clump masses. It 
however, requires that the instantaneous SFE, i.e. the mass of 
stars formed from a particular mass of de nse gas, acting on eac h 
clump varies as a func tion of clump mass dUrquhart et alj|^14bt) . 
dUrquhart et alj[^14bh suggest that the SFE is lower for the high- 
mass and high-luminosity sources in which the high-mass mem¬ 
bers of a cluster form first. The disp arity between the obse rved 
and cluster relationship s is seen by both lSridharan et aDd2002h and 
lUrauhart et alj d2014bh . Both of those studies show that the lower 


mass clumps are over luminous, whilst the higher mass clumps are 
not as luminous as one would expect if each source was expected 
to be the site of cluster formation and a standard IME was sampled 
at a constant SFE. However, the observed dependence is consistent 
with the relationship between accretion luminosity and stellar mass 
{L oc and that of accretion luminosity and final stellar mass 

(LocM^'^°: lMcKee & Tanll2003h . This implies that these sources 
are either still accreting or have just reached the final stellar mass. 


6 SUMMARY & CONCLUSIONS 

We have combined a series of surveys tracing infrared YSOs with 
the Bolocam Galactic Plane Survey (BGPS) sample of star-forming 
clumps within two slices of the Galactic plane, in order to inves¬ 
tigate how the star-formation efficiency, measured as the ratio of 
infrared luminosity to clump mass, varies with Galactic environ¬ 
ment. The two slices of the Galactic plane, centred at Galactic lon¬ 
gitudes of / = 30° and / = 40°, cover spiral-arm regions, interarm 
regions and the end of the Galactic Long Bar, which contains the 
star-forming region, W43. 

As a result of positionally match ing infrared sources to B GPS 
clumps using the sample subsets from lEden et al.l ( |2012L|201^ . we 
find a total of 348 candidate YSOs associated with 298 distinct 
clumps, which accounts for 29 per cent of the BGPS sources found 
in the two regions. The distances, and Galactic environments, of the 
BGPS sources and the luminosities of t he associated YSO sourc es 
were taken from the derived distances in lEden et alj ( 120121120131 . 

In order to obtain a YSO luminosity independent of the clump 
mass, we used only wavelengths shorter than 70 qm to calculate 
Lir and the 1.1-mm flux density for Mdump - The method for cal¬ 
culating the luminosities was to use a trapezium rules integration 
in that wavelength range as the use of the bolometric luminosity 
from the Robitaille SED fitting tools dRobitaille et alJl2006Ll200'T) 
included the contribution from the envelope in which the YSO is 
embedded, i.e. the clump, and an independent measure of the YSO 
luminosity from the clump component was required. The absolute 
luminosities found using both methods were found to be different 
but the relative differences were consistent across all distances and 
environments. 

We compared the ratio of infrared luminosity, Lir, to the 
clump mass, Mdump, as a function of Galactic environment, and 
the properties of clumps hosting a candidate YSO and our main 
findings are: 

(i) In the 1 = 30° region, the Lin/Mdump ratio is larger in the 
background spiral arms by a factor of ~ 2.5 compared to the 
other / = 30° populations, including the W43 star-forming re- 
gion. This is consistent with the results of iMoore et al.l ( l2012h 
who found that the Sagittarius and Perseus arms had a higher 
Lin/Mdoud ratio than the Scutum-Centaurus arm. This may 
imply that the future SER of W43 wi ll be higher than the 
curren t observed rate, as suggested by iNguven Luong et alj 

( l201lh . 

(ii) The luminosity distributions for the Scutum-Centaurus tan¬ 
gent region and a combined spiral-arm sample, made up of 
Sagittarius and Perseus arm sources, are consisten t with each 
other. This supports the Mdump/Mdoud results of lEden et alj 
|2 oI 3) that the mass function at the end of the bar is not sig¬ 
nificantly different from any other spiral arm in the Galaxy 
and that there is no enhancement in the SFR per clump mass 
due to the interaction of the bar with spiral arms. 
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(iii) The observed distribution of -LiR/Mdump ratios for all 
clumps is consistent with a log-normal probability function, 
implying that those with extreme values, both low and high, 
are the result of the same combination of random processes 
and do not indicate the presence of fundamentally different 
conditions. 

(iv) Assuming a single star forming from each YSO, the average 
absolute SFE, ratio of stellar mass to clump mass, was found 
to be 1.9 ± 0.2 per cent. 

(v) The ranges of L/M found for the individual clumps, 

molecular clouds and spiral arms were found to 
be 0.04 ±0.01-140 ±22, 0.02 ± 0.01-68 ± 13 and 

1.00 ± 0.02 - 2.9 ± 0.2 Lq/M© respectively. 

(vi) Combining these results with others investigating the role of 
spiral arms on the star-formation process it is found that, on 
kiloparsec scales, there is little evidence of significant varia¬ 
tions in SFE that can be attributed to the action of the arms on 
the ISM. Differences in average Lm/Mciump found in spiral- 
arm sources may be due to gradients in the SER or local vari¬ 
ance in the limited samples. The most significant variations 
in SEE are found to be on the scales of individual molecular 
clouds and clumps. 

(vii) There is a linear relationship between the Lir and 
Mciump with a best-fit power-la w exponent of 0 .99 ± 0.04, 
which is consistent w ith that of iMolinari et alj ( l2008h and 
lUrauhart et al.U2014bl) . 
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